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1 Introduction 



Non-leptonic weak decays of B mesons are very interesting for several reasons. First, CP 
violation in the B-meson system will eventually give us information about the CP violating 
phase in the Cabibbo -Kobayashi -Maskawa (CKM) mixing matrix JIJ]. Second, non-leptonic 
weak decays will give additional clues for determining the absolute values of the quark 
mixing parameters, in particular the ratio \V u b/V c b\, although it is expected that more solid 
information will come from semi-leptonic B decays. Last, the dynamics of the non-leptonic 
weak decays in the framework of the standard model is only poorly understood. One of the 
problems in calculating the transition amplitudes for non-leptonic weak decays is that one 
needs to evaluate the hadronic matrix elements of certain four-quark operators which can 
be done in QCD only with non-perturbative methods. The usual route to calculating these 
hadronic matrix elements is to start from the effective, QCD corrected, Hamiltonian for the 
Ab = 1 non-leptonic decays in the six-quark model, i.e. including the t quark || ||], |J. This 
gives the weak Hamiltonian in terms of four-quark operators. For computing the hadronic 
matrix elements of these four-quark operators the factorization approximation is used 

[7). Then the hadronic matrix elements are given in terms of current matrix elements 
(matrix elements of two-quark operators) as they appear also in semi-leptonic decays. These 
current matrix elements are much easier to calculate and many models || |9|, have been 
proposed for them. When these current matrix elements are approximated by a single pole 
based on the idea of vector dominance we arrive at the pole- dominance model of two-body 
non-leptonic decays. This model has been applied to the calculation of non-leptonic decays 
of charmed mesons JTTJ. Many years ago, this model was used already for the discussion of 
strange particle decays. In particular, it is known, that this model provides a basis for a 
description of K — > 2n decays |T2| . 

In this work we apply the pole- dominance model to two-body non-leptonic B decays and 
show its strong relationship to the usual factorization model |7], Therefore this latter 
model will be our starting point where we also use as input the parameters of the current 
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matrix elements as given in 0. This allows us to extract various strong coupling constants 
between B mesons, their vector and axial- vector counterparts, and light and charmed mesons. 
These coupling constants and also the form factors of the current matrix elements will be 
compared with other information coming from QCD sum rule and QCD lattice calculations. 
For simplicity, we consider only the lowest lying poles. The inclusions of higher lying poles 
would make the extraction of the strong coupling constants impossible. 

The outline is as follows: In section 2, we explain the relationship of the usual ap- 
proach with the pole-dominance model using as example a collection of B decays into two 
pseudoscalar-mesons (B — > PP). We calculate the two-body branching ratios in parallel 
for the usual approach and two versions of the pole-dominance framework. In one version 
we calculate the residues of the poles on the mass shell of the intermediate state. This 
reproduces the usual approach. The other version follows essentially the work of Bedaque 
et al.|ll|| for D decays where a particular off-shell extrapolation has been adopted. These 
results are compared to experimental data when they are available. In section 3, we proceed 
to the more complicated B — > VP and B — > VV channels, where V are light or charmed 
vector mesons. Section 4 is reserved for a summary and some conclusions. 



2 B decays into two pseudoscalar mesons 

To start, we choose as an example for the decays B — > PP, the channel B° — > 7r + 7r~. The 
relevant effective weak Hamiltonian can be parameterized as 

n eff = ^ + c 2 o 2 )] + h.c, (i) 

where the Cj (i=l,2) are the short-distance Wilson coefficients defined at a scale /i of the order 
of the heavy quark mass m?, and the Oi are the local quark operators with the appropriate 
quantum numbers. For simplicity we have neglected the strong and electroweak penguins. 
The electroweak penguins are found to give small contributions to the mode B° — > 7r + 7r~. 
The strong penguins could be included easily |T4 . 
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The explicit expressions of the operators 0\ and O2 are: 

O x = (d a ^Lup)(u p j^Lb a ), (2) 
2 = {da^Luajiu^^Lbp), 

where a,j3 are colour indices and L = |(1 — 75). For the Wilson coefficients at the scale 
H = m b = 4.8GeV we use the values ||14| : 



d = -0.324 C 2 = 1.150. (3) 
These coefficients are regularization scheme independent at next-to-leading logarithmic 



precision as obtained by Buras et al.|15| for = 350MeV. With the factorization hypoth- 



esis and after Fierz reordering we obtain from ([!]) and (0) 

< n+n-\H eff \B° >= ^=-V ub V* d (±d + C 2 ) < 7i+\u lfl Lb\B° >< 7r-\d^Lu\0 >, (4) 

where N is the number of colours. 

The current matrix elements in (^) are evaluated in terms of form factors F\ and Fq for 
the decay B° — > and the pion decay constant /„- = 0.132GeV as follows: 

,2 _ ^2 
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<n + \u^Lb\B° >= l - 



m o - m_ 
(Pb+P*)» ~ 2 % 



1 2 2 

1 m B — mi 

2 ~f 



<Tr-\d^Lu\0>=~Uq» (6) 



where q = Pb — Pn- With these definitions we get the following result for 



< n+n-\H eff \B° >= i^V ub V^ (Id + C 2 ) U(m% - ml)F B ^(ml) (7) 

In (||) Fi(Fq) are the form factors of 1~(0 + ) (transverse (longitudinal)) projections of the 
vector current. At q 2 = 0, we have -F (0) = -Fi(O) to cancel the pole at q 2 = in (|5]). But Fi 
and F differ for arbitrary q 2 . In ([?]) we need F^^ n only for q 2 = m 2 , which is small, so that 
we can set Fq^ u = Ff^ n . Of course for F^^ 7T (q 2 ) = Ff^ n (q 2 ) with arbitrary q 2 in (||) we 
obtain the same result. The 1 _ (0 + ) part of the current (u^f^Lb) has the quantum number 



of the B*~ (Bq) meson with mass mg. = 5.32GeV (m# = 5.73GeV) (We denote the 3 Pj 
bound states by Bj). If according to the vector-meson-dominance model (VDM) the form 
factor F^ 71 is dominated by just this 1~ resonance with mass mg. we have 

m B* ~ Q 

In (^j) the current coupling f B * is defined by the vector current matrix element 

< 0|M7 M 6|fi* >= fB-rriB*^, (9) 
and the strong coupling is defined through the Lagrangian 

C-B'B. = gB*B,B;B &> 7T. (10) 

Similarly we define the current matrix element of the scalar resonance state Bq\ 

<0\uYb\B >=if Bo P Bo , (11) 
and the coupling of Bq to B n + by 

1 ^ 

£-B Bn = 9B Bnd^B B <9 M IT, (12) 

9B Bn is dimensionless. With these definitions the form factor F^^ n (q 2 ) is given in the 
scalar-meson-dominance model by the formula 

F ° {q)= ml -q* ■ (13) 
With ([13|) the weak transition matrix element (|7D can be written as 

< n+n-\H eff \B° >= i%V uh V: d ()-C x + C 2 ) U ™} ~ m \ g Bo B.m B J Bo . (14) 

V I \iv / m Bo — m n 

Here 9b q b-k and /b are on-shell couplings to the scalar resonance Bq which is on the mass 
shell. This means that the off-shell extrapolation factor of the current matrix element is 
given by just a single pole in flT3"|) . 

The result ([14]) for the transition amplitude B° — ► 7r + 7r _ can be obtained more directly 
from the pole diagram in Fig. la. In this diagram, the initial B° can go into a B*~tt + or 
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(a) (b) 
Figure 1: Feynman diagrams for B° — > 7r + 7r~ in pole model. 

Bq7t + via the strong interaction couplings ( |TD| ) or (|T2|). The or B converts into a 
7r~ via the weak interaction (]!]) shown in terms of quark lines in Fig. lb. The weak matrix 
element is evaluated in the vacuum insertion approximation with the following result for the 
intermediate Bq\ 

< 7r-\n eff \B- > = ^V ub V: d (^d + C a ) < n-\dYLu\0 >< 0\u llM Lb\B- > 

= -^V ub V: d ^C 1 + C^f^m% o f Bo . (15) 

In (|15D p n ~ = p B - and p 2 B - = vn 2 ^ and the current matrix elements in (|6|) and (|TT|) have 



been used. Inserting the propagator of the intermediate B meson and the strong vertex 
the result for the decay matrix element is: 



< n+n-\H eff \B° >= i%V ub V: d f^d + C 2 ) f^T^HLg m f (i 6 ) 

V2 \N J m Bo - 7774 

The result (H) agrees with (|1J) as to be expected. To achieve this, two ingredients were 



essential. First, the coupling Bq — > Bit in ( |T!?| ) has such a form that it vanishes for m? B = m 2 , 
i.e. the coupling of a scalar meson to two pseudoscalar mesons of equal masses would be 
different. Second, in (113) we have set p n - ■ p R - = p 2 = m\ 1 . This means that the weak 
transition matrix element Bq — > 7r~ is evaluated at the on-shell point of the intermediate 
Bq resonance, whereas the propagator is calculated at (p B o — p n +) 2 = v\- — m 'i- This is in 



the spirit of the meson-dominance approximation and agrees with the result (p~3|) where the 
form factor F B ^ 7r (g 2 ) is written down directly in this approximation. 



In the pole model calculations of ref. |TT[ a different route has been taken. Instead of the 
Bq state the intermediate state is the B* with spin 1~. In this case the pole cancels by the 
factor m 2 B , — q 2 in the numerator and the result is 

< n+n-\H eff \B° >= i%V ub V: d + C 2 ) /„ (m * " ^ g B * B7r f B * ■ (17) 



The difference between ([[D and fll7|) is small since Ff^(0) = F^ n (0) ~ F^K). It 



is clear however, that in this evaluation the residue of the pole l/(m B * — q 2 ) vanishes like 
(m B * — q 2 ), so that the pole is canceled. This is apparent also from ( |l5f) when we evaluate 
< n~\H e ff\B*~ > in the vacuum saturation approximation. Then we obtain the factor 
e • p n - = e ■ p B * where e is the polarization vector of the B*. Of course, this factor vanishes 
on the mass shell. We think that our evaluation of the pole model is more appropriate since 
the pole is kept, so that our result is correct for q 2 — > m Bo . Since a large extrapolation from 
q 2 = m 2 Bo to q 2 = m 2 is involved, one must be careful in the calculation of the numerator. 
From fljjD we have 

_ F o B -(0)m Bo 

9b q Btt — t • { L ° ) 

JBq 

and a similar equation for g B * Bn following from (g). 

Ff~* n (0) = Fq^ 7 "^) has been calculated in various approaches: quark model, QCD sum 
rule and QCD lattice calculations. Thus either from these results or in case the branching 
ratio for the decay B° — > 7t + tt~ becomes precisely known experimentally the relation (16) 



allows us to calculate the coupling constant g BoBn , if f Bo is determined from other sources. 
Similarly g B * B7T is calculated from Ff^ n (0), where we need in addition the current coupling 
f B *. Since the experimental branching ratio of B° — > 7r + 7r _ is not very well known yet (see 
Table 5) we must rely on F 1 B ^ 7r (0) = F^ _>,r (0) as calculated in various models. We shall do 
this later and compare to direct theoretical results for g BoBlT or g B * B7T . 

Our result for the weak transition matrix elements has the factor a 2 = Ci + C\/N, which 
should be evaluated with N=3. In the following we shall evaluate them with N=2 giving 
a 2 = 0.988 using (|J). a 2 depends on N only mildly, since |Ci| < C 2 , so that choosing N=3 
would give similar results for all colour unsuppressed decays, for which the transition matrix 



elements are proportional to 02. In the following we shall denote this class of decays by class I. 
The class II stands for the colour suppressed decays which have matrix elements proportional 
to ai = C\ + C 2 /N, which depends strongly on N. The third class is a superposition of class 
I and class II matrix elements. Recent comparisons of experimental data for many two-body 
B decays show that on the basis of the factorization hypothesis and BSW matrix elements 
from |13|] a reasonable fit to these data gives N ~ 2 fTEfl . N=3 would produce very bad 



results for the measured class II decays and the class III decays can be fitted only with 
a positive a±. For example, for N=2, 3 or 00 we have the following values for a x using 
(g): ai = 0.251,0.059,-0.324. The choice N = 00 is favored in D decays 0. The fact, 
that the choice N=3 is not possible points into the direction that so-called non-factorization 
contributions are significant, in particular in the class II and class III transitions. Of course, 
it is not certain, that these contributions can universally accounted for by choosing N=2. It 
is more likely that these contributions depend on the particular decay channel. For fits and 
constraints on the non-factorization terms in various channels see [T7). We also emphasize 
that the result iV ~ 2 from the fit in [[16] depends on the BSW |13[ choice of transition 
form factors and current coupling constants. We shall assume N = 2 independent in which 
pole model the weak matrix elements are evaluated. This can be justified empirically since 
the magnitude and also the sign of ai is deduced from such final states, where the two pole 
models coincide, as will be seen later in our presentation for the PV and VV results. If one 
gives up the requirement that the should be independent of the regularization scheme one 
can find short distance coefficients C\ and Ci in next-to-leading logarithmic precision which 
give ai = 0.2 with N=3 by selecting the appropriate scheme [|H|. 

Other class I decays going into two pseudoscalar mesons are B° — > D~D + , D~D + , 
K~D + , 7i~D + and B~ — > D~D°, D~D°. These decays are calculated from H e //'s where 
u — > c and/or d — > s with the same short distance coefficients as in (3). Of special interest 
is the decay B° — > tt~ D + which belongs also to the colour unsuppressed class I. This decay 
can proceed also through the annihilation diagram. The diagrams in the usual factorization 
approximation for the quark diagram and for the pole model are shown in FigJ^a and ^0. 
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Figure 2: The diagrams of B° — > D + n in the usual factorization approximation for the 
quark diagram (a) and for the pole model (b). 



In the pole model according to ref. [11], the transition matrix element corresponding to 
the diagram in Fig.[|b is 

< n-D + \H eff \B° >= - % %V ch V: d (±C 2 + C x ) /bo (m ^7 m * ) g D .o D+w -f D .m D .. (19) 



V2 



With M) this can be written as 



< TT-D + \n eff \B° > 



N 



G p 

' : 7f 



rn 



D* 



V cb V: dai f B0 (m 2 D - ml)F?^{U). 



(20) 



Compared to the dominant pole contribution, which is given by the analogous equation to 



2 

fll4|) , the annihilation contribution is suppressed by the factor ^r-^ which is approximately 
0.03. 



The fact that in ( 2D ) the pole l/(q 2 —m 2 D ,) has canceled against the same factor (q 2 —m 2 - ) *) 
in the numerator indicates that the residue of the pole vanishes on the mass shell of the D*° 
resonance. This is already obvious when one evaluates the weak transition matrix element 
in the vacuum saturation approximation 



< D*°\H e ff\B° >-- 



AG F 
V2 



V cbV* d a x < L>*°|^y M Lu|0 >< 0\d^Lb\B° > . (21) 



Table 1: Values of pole masses in GeV. 



Current 


m(O-) 


m(l-) 


m(l+) 


m(0+) 


dc 


1.8693 


2.010 


2.4222 




uc 


1.8645 


2.0067 


2.4222 




sc 


1.9685 


2.1124 


2.53535 




ub 


5.2789 


5.3248 


5.37 


5.73 


db 


5.2792 


5.3248 


5.37 


5.73 


sb 


5.3693 


5.41 


5.9 


5.89 


cb 


6.264 


6.337 


6.730 


6.700 


cc 


2.9798 


3.09688 


3.51053 


3.4151 



The first current matrix element in fl2~l"|) is proportional to the polarization vector of the D*° 
whereas the second current matrix element is proportional to p B o = Pp*o, so that the product 
vanishes. This indicates, that the pole-model calculation of the annihilation contribution is 
also ambiguous. 

If we start from the quark diagram in Fig.||]a and do the vacuum insertion approximation, 
but then evaluate the matrix element < 7r~ D + \c^f fl Lu\0 > directly in terms of the form factor 
F ( f _ * 7r (g 2 ) we obtain the usual result: 

<n-D+\H eff \B°> = -&V a V: d aJ B0 (m 2 D -ml)F D ^(m 2 B ) (22) 



V2 m 2 Do - m% 



This means that the annihilation contribution is suppressed further by the additional factor 
m D /( m B ~ m D )i "which results from the fact that in (^) the current transition form factor 
is highly off-shell at q 2 = m B , whereas the off-shell extrapolation in the pole diagram is 
determined by the mass of the D*°. It is clear that the pole-model calculation starting from 
FigJ^a with a + intermediate state leads to the same result (|22|). 

Such an annihilation contribution is in principle also possible for B° — > 7r + 7r~. But it 
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vanishes since the particles in the final state have equal masses. 

We conclude that the pole-model calculation of the non-annihilation diagrams as depicted 
in Fig. la agrees with the usual calculations of the B — > PP decays where the form factors are 
approximated by a single pole (compare ([TJ) with (|T^) ). If a different off-shell extrapolation 
as advocated in ref . [11[] is assumed we obtain the result ([H]) which differs actually very little 
from (|I6D since F^ w (0) = Ff^ w (0) and the extrapolation of F^^ n (q 2 ) from q 2 = to q 2 = 
m 2 has only a small effect. Of course, this is different for decays like B° — > D~D + , D~D + , 
for which the form factor F^~ >D (q 2 ) must be computed for q 2 = m 2 Dl which differs from 
Ff^ D (0) to be inserted into flT7|). 

Class II decays are B° — > 7i tt°, tt°D°, K°r] c , and B~ — > K~7] c . The decays B~ — > 
7t~7t°, tt~D° and K~D° are class III decays, where contributions of colour suppressed and 
unsuppressed matrix elements interfere. 

To obtain an overview on expected branching ratios for all these decays we shall present 
results in two schemes: (i) the pole-model calculation in the form (|T4| ) which is equal to the 
usual quark diagram computation in terms of form factors F (0); (ii) the pole model result 
flTTD where the pole canceled and the result depends on -Fi(O). In scheme (i) the annihilation 
contribution which occurs in some channels is negligible, either due to strong suppression in 
the quark diagram result fl22|) or due to the same suppression for the pole model. In scheme 
(ii) we evaluate the annihilation diagram with the formula (|20|) which involves a different 
off-shell extrapolation of the residue of the pole. It is clear that in this form the annihilation 
terms are larger. The scheme (i) will be denoted BSW model since the pole model calculation 
agrees with the usual computation and the scheme (ii) will be denoted pole model. 

For these calculations we need the values of the masses of the particles in the initial and 
final states and of the pole positions together with the pseudoscalar (/m) and vector meson 
(fv) decay constants. The masses are taken from the PDG tables |19[ or from the table in 
[0. The decay constants fw and fy for pseudoscalar and vector mesons, respectively are 
taken from the following sources. /„-, fx, f P and fx* are taken from 0. fo s = 288MeV 
corresponds to the CLEO measurement |2(J of the leptonic decay of the D~ . We assume 
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Table 2: Values of decay constants. 



Particle 


fu (MeV) 


Particle 


fv (MeV) 


7T 


133 


P 


221 


K 


162 


K* 


221 


D 


250 


D* 


250 


D s 


288 


D* s 


288 


B 


140 


B* 


160 


B s 


158 


B* s 


181 


B c 


360 


b; 


320 


Vc 


393 


j/4> 


382 



Id* = fD s and fix by the ratio foj /d = 1-15, which is the average value of this ratio 
obtained in various lattice calculations [pTJ] . fjj = fr>* is also fixed by assumption. f B and 
]b* are from QCD sum rule calculations in [p2| . The values for fg a and are fixed by the 



ratio Jbs/Jb = fs*/fB* = 1-13 from pl| . The values for fg c and fs* have been found in |23| . 
The ratio of f^J fj/^ = 1.06 is found in [^4] and fj/^ is determined from the leptonic width of 



the J/ip as in Of course, most of these values for the decay constants are not definite yet 
and somewhat other choices can be made. This information is collected in Table 1 (masses) 
and Table 2 (decay constants). The form factors F\(0) = Fq(0) are written down in Table 3, 
also taken from [[/]]. The CKM matrix elements are given in Table 4, together with the B° 
and B~ lifetime data [|19j. Our results for the branching ratios in the pole model and the 
BSW model are collected in Table 5. In this table we also give the amplitudes for the BSW 
and for the pole model which show in particular the contributions proportional to a\ and a2 
for the class III decays. For some of the channels the branching ratios have been measured. 
The world-average values are reported in p^| . For other channels some new measurements 



and upper limits exist f25|. These data are given in the last column of Table 5. 



From the results in Table 5, we observe the following. Pole model and BSW results lie 
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Table 3: Form factors at zero momentum transfer . 



Decay 


TP J? 

Pi — Pq 


V 








D - 


-* 7T 


U.oyz 










u - 


■* K 


U. ( oz 










u - 


P 




1 £ 

l.zzo 


U. M O 


U.yzo 


U.ooy 


D* - 


— > 71 






n 77 £ 
U. ( / o 


U.yzo 


u.ooy 


D s - 


-> if 


0.643 










B - 


-> 7T 


0.333 










B - 


-> if 


0.379 










B - 


-> D 


0.690 










B - 


^ P 




0.329 


0.283 


0.283 


0.281 


B - 


► if* 




0.369 


0.328 


0.331 


0.321 


B - 


> D* 




0.705 


0.651 


0.686 


0.623 



Table 4: CKM matrix elements and life time of B mesons. 



Name 


\Vud\ Wcs\ Kis Kb \Vub\ 


T B o (s) 


t b- (s) 


Values 


0.975 0.9742 0.221 0.221 0.038 0.0035 


1.56 x 10" 12 


1.62 x 10" 12 
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Table 5: B — ► PP Amplitudes(10 7 GeV) and Branching Ratios in two versions of pole 



models and experimental measured branching ratios or upper limits (90%C.L.) [16 







Pole model 




BSW model 




Experimental 


Channel 


Amplitude 


BR 




Amplitude 


BR 




BR 


B- - 


-» D-D° 


2.91 a 2 


5.44 x 10" 


-4 


3.16 a 2 


6.39 x 10" 


-4 




B~ - 


- D;& 


14.8 a 2 


1.36 x 10' 


_2 


16.2 a 2 


1.63 x 10" 


_2 


(1.36 ± 0.43)10~ 2 


B°- 


-> 7r~7r + 


0.35 a 2 


1.05 x 10' 


O 


0.35 a 2 


1.05 x 10" 


o 


(0.7 ± 0.4)10~ 5 


B°- 




1.89 a 2 


2.68 x 10" 


-4 


1.90 a 2 


2.71 x 10" 


-4 




B°- 


-> 7T-P + 


6.84d2-1.03oi 


3.29 x 10" 


-3 


6.84 a 2 


3.56 x 10" 


-3 


( 3.1 ±0.4)10~ 3 


B°- 




14.8 a 2 


1.31 x 10" 


-2 


16.2 a 2 


1.57 x 10" 


-2 


(0.74±0.28)10~ 2 


B°- 


-> P~P + 


2.91 a 2 


5.22 x 10" 


-4 


3.16 a 2 


6.14 x 10" 


-4 


< 1.2 x 10~ 3 


B- - 




12.6 oi 


6.13 x 10- 


-4 


16.9 a x 


1.11 x 10" 


-3 




B°- 




12.6 ai 


5.90 x 10" 


-4 


16.9 oi 


1.07 x 10" 


-3 




B°- 


+ 7rV 


0.17 ai 


3.39 x 10" 


-7 


0.17 ai 


3.39 x 10" 


-7 


< 9.1 x 10~ 6 


B°- 


+ 7T £> 


5.01ai-0.72ai 


9.04 x 10' 


-5 


5.60 a x 


1.54 x 10" 


-4 


< 4.8 x 10" 4 


B- - 




1.89a 2 +1.81ai 


4.32 x 10' 


-4 


1.90a 2 +2.02ai 


4.55 x 10" 


-4 




B- - 




6.84a 2 +7.08ai 


5.91 x 10' 


-3 


6.84a 2 +7.92ai 


6.20 x 10" 


-3 


(5.0±0.5)10~ 3 


B- - 




0.25a 2 +0.25ai 


8.57 x 10" 


-6 


0.25a 2 +0.25ai 


8.58 x 10" 


-6 


(9±|) x 10" 6 
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Table 6: The form factor Ff^ n 



0) at q 2 = in different models. 



Reference method 





in 





29 



3C 



31 



32j 



in 



34 



a 
a 
a 
b 
c 
c 
c 
c 
c 
c 

d 

c 
c 
c 



Ff^(0) 



0.333 
0.293 
0.21 ±0.02 
0.53 (0.89) 
0.29±0.01 
0.26 ±0.02 
0.26 ±0.02 
0.23 ±0.02 
0.4 ±0.1 
0.24 
0.37 ±0.05 
0.10-0.49 
0.23-0.43 
0.21-0.27 



a Quark model 

b HQET and CPTT 

c QCD sum rules 

d experimental value from B(B° — > tt~1 + i>) 
e Lattice calculations. 
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very near together. Exceptions are the decays B~ — > D~D°, D~D° and 5° — > D~D + , 
D~D + which differ due to the extrapolation of F^~^ D from g 2 = to q 2 = m 2 D . The 
annihilation contribution in B° — > n~ D + has little effect. However, it is significant for the 
decay B° — > tt°D , where it is proportional to the same Wilson coefficient a\ as the direct 
contribution. Except for the decay B° — ► D~D + we observe very good agreement between 
the calculated and the measured branching ratios. Both pole model versions agree with the 



data, so that they can not be utilized to exclude for example the model based on ref.[fLl . 

For the results in Table 5 we needed the values of three form factors F^ D , and 
Ff^ K taken from Table 3. Since only for four channels, namely B~ — > D~D°, ir~D° and 
B° — > D~D + , 7i~D + , experimental branching ratios exist, we gain only information on 
Ff^ D (0) and Ff~* n (0). Since our results agree quite well with the measured data we are 
confident that the BSW values for these form factors are reasonable. In this connection it is of 
interest to search in the literature for other theoretical predictions on these form factors. Such 
predictions based on various methods are collected in Table 6 for F 1 B ^ 7r (0) and in Table 7 
for Ff^ D (0). The values for Ff^ 77 with two exceptions are very similar. Calculations using 
similar methods give approximately the same result. Furthermore the values predicted from 
QCD sum rules and constituent quark models of various kinds are very near together and 
agree with the average value obtained from QCD lattice calculations. From the comparison 
in Table 6 we are confident that Ff^ n (0) ~ 0.3 is a reasonable value which agrees also with 



the recent measurement of the branching ratio for B — > 7T l + v [34]. Fj ^(0) = 0.37 ±0.05 



is deduced from the measured branching ratio, using \V U b\ = 0.0035 from Table 4 and the pole 
masses from Table 1. This value depends on the way the q 2 dependence of the form factor is 
parameterized. We assumed for consistency the approximation with one pole. Results with 
other form factor assumptions can be found, for example, in pUf . 

The situation is similar for the form factor Ff~* D (0). The results obtained from quark 
models agree quite well and agree nicely with the value calculated from the measured branch- 
ing ratio B(B° — > D~l + v{) = (1.9 ± 0.5) x 10 -2 fl9 |. Of course, the experimental value 



F^^ D (0) = 0.71 ±0.10 and also the results obtained from the Isgur -Wise normalization |39| 
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Table 7: The form factor Ff^ D (0) at q 2 = in different models. 



Reference method 



H 





n 



a 
a 
a 
c 
f 



0.690 
0.684 
0.63 
0.62 ±0.06 

0.74 
0.71 ±0.10 



f HQET with monopole form factor and £(1) = 0.91 [|0 
g experimental value from B(B° — > D~l + v). 



in the HQET depend on the assumed q 2 variation of F^~* D . The results in Table 7 under 
f and g are obtained with single-pole approximation with masses from Table 1. Somewhat 
smaller values as from the extrapolation with the single-pole formula are obtained with the 
usual exponential extrapolation to q 2 = on the basis of slope calculations from QCD sum 
rules or lattice calculations or from fits to the experimental data of B — > D*lu decays p0"| . 
Taking into account that the Isgur-Wise normalization is £(1) = 0.91 ±0.04 |4(J due to mass 
and perturbative QCD effects one arrives at Ff^ D (0) ~ 0.61, instead of F^~* D (0) = 0.71 in 
Table 7. We conclude that the BSW value is a good average. 



We are now in the position to calculate gB*Bir, 9b*bk and Qb*bd from flTq). To be 
definite we use the BSW values of the corresponding form factors from Table 3 together 
with the decay constants from Table 2. We obtain for the strong VPP coupling constants 
9b*b-k = 11.1, 9b*bk = 11-7, and Qb*bd = 13.7. Surprisingly, these three couplings lie 
close together, so that the spread of the form factors Ff^ p (0) (P = 7T, K, D) is mostly 
due to the different values of fs*, Jb* and fs* and the different mass values. If we take 



fs* = 0.400GeV, i.e. the upper limit of the QCD sum rule result p3[ we obtain Qb*bd = 11.0. 



In potential model calculations fs* comes out even larger: fs* — 0.5GeV [|4|. Results for 
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9b* b-k based on QCD sum rule calculations or on the quark model with chiral HQET are 
collected in p2]] . Our result g B * BlT = H I can be compared with various QCD sum rule 
calculations. The calculations use different approaches. They are either based on expansions 
near the light-cone or use the short-distance expansion in connection with the soft pion 



limit. Results obtained with the first method are gB*Bw = 14 ± 3 P^fl, g B * BlT = 10 ± 2 [12 



Qb'B-k = 10.4 ± 2.0 |4"3"| , and g B * BlT = 7± 2 J44|]. The second method has given the following 
results: Qb*b-k = 14.5 ± 1.5 J^] an d Qb*Btt = 11-2 ± 2.0 [[4^]. These results depend on the 
values obtained for f B and f B *, which are also taken from QCD sum rule calculations. In 
[||, H, (H these constants are f B = (150±20)MeV and f B * = (190±10)MeV. This explains, 
for example, the different results in p2j and with the light-cone sum rule. In the original 
result for f b fB* 9b* b-k-, as it follows from the sum rule, the two evaluations gave the same 
result, namely fsfs'SB'Sir = 7.9 x 10~ 3 GeV 2 . The smaller value of gB*Bn, obtained by 
Narison et al. |44]] is due to the fact, that these authors included 0(a s ) corrections to the 
heavy meson decay constants which result in larger values for f B and f B * ■ Taking all these 
uncertainties into account the value of g B * Bn obtained from the non-leptonic decays agrees 
reasonably well with the QCD sum rule results. 

To our knowledge there exist no QCD sum rule calculations for g B * B R and g B * BB to 
which we could compare our results. 

Actually in the BSW pole model the decay amplitude depends on the form factor Fq ^(0) 
(see ©) not Ff^ n (0). At q 2 = we have F B_>7r (0) = Ff^(0) which yields an equation for 
g BoB7T in terms of g B * BlT : 

m B fs* , oq n 

g bqBti — -j. — gB*Bn y^o) 

m B* /b 

i- e - g BqBtt — gB*Bir- Inserting masses and assuming f Bo = f B , we obtain g BoBn = 13.7. This 
result can be compared with a recent QCD sum rule calculation of this coupling constant 
H|. Adjusting their result to f B = f Bo = 0.14GeV and to our definition of the BqBtt 
coupling in ( |12|) the result in |46| is gB B-w = 4.5 and gB Q B-w = 6.3 depending on the QCD 
sum rule method used. These values are more than a factor of two smaller than our result. 
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Figure 3: Feynman diagrams for B° — > tt + p (a) and for B° — > 7r p + (b) in pole model. 

It is clear that these sum rule results are not compatible with the relation (|23"D . Since our 
coupling (|l2|) and the BqBtc coupling used in have different off-shell extrapolations it is 
also not clear to us whether this has any effect on the sum rule evaluation. 

We conclude that the coupling constant gB*Bn obtained from non-leptonic decay data 
is in reasonable agreement with QCD sum rule results. The coupling constants Qb^bk and 
9b*bd have rather similar values to gB*Bn- It would be interesting to know whether these 
relations can be explained in the framework of QCD sum rule calculations. 

3 B decays into pseudoscalar and vector and two vec- 
tor mesons 

First, we consider decays of B mesons into one vector and one pseudoscalar meson: B — > VP. 
As examples we choose B° — > 7i + p~ and B° — > 7i~p + , the analogous channels to the decay 
considered in the previous section. The pole model diagram for B° — > tt + p~ is shown in 
FigJ^a. The intermediate state is the B* resonance which decays weakly into a p~~ meson. 
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The result of the pole model evaluation is 



< vr> \Heff\B >= — =V ub V ud a 2 m p f p e ■ (p B + p n ) = — , (24) 



which can be written as 

< n + p-\H eff \B° >= ^V ub V: d a 2 m p f p e* ■ (p B + p„)F?^ (m 2 p ), (25) 

using the vector-dominance relation (^) for the form factor F^ 71 . This result is identical to 
the result of the usual quark model calculation of BSW. It is clear that an intermediate scalar 
resonance Bo gives a vanishing transition matrix element. So, for this case the two pole model 
calculations yield a unique result. The same formula applies to the decays B° — > D + D*~, 
D+D*-, D+p- and D+K*~. 

For the decay B° — > 7r~p + , we need in the usual calculation from the quark diagram 
the axial- vector current matrix element between B° and p + . This is parameterized by three 
invariant form factors A , A\ and A 2 [§: 

< p + |m7m75&| B° > = z fe* - e -^-qA (m B + m p )A^ ~* p (q 2 ) 



— — e • q) — - — 

q z J m B + m p 



i\(pB + v P )ii - 2 — -q» (e* • q) 



+ i^!^), (26) 
q z 

where q = p B — p p = p n and e* is the polarization vector of the p + . To cancel the poles at 
q 2 = 0, we must have 

2mX" P (0) = (m B + m p )Af-* fi (p) - (m B - m p )A^ p (0). (27) 

With (|26|) we obtain for the weak transition matrix element the usual result: 

< n-p + \H eff \B° >= ^£v ub V^a 2 f n (e* ■ p n -)2m p A^ p {ml). (28) 

The form factor Aq^ p can be approximated by a single pole with spin 0~. Analogously to 
(j§D and (^) this meson-dominance approximation yields 

A B ^P(n 2 \ - 9BBpfBm 2 B 

° iq) ~ (m 2 B-q 2 )m P { ) 
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The intermediate state is the well known B~ particle with mass m B . When we evaluate the 
pole-model diagram in Fig.^D with an intermediate 0~ particle directly we obtain the same 
result, namely 

< n-p+\H eff \B° >= ^V ub V: d aMe* ■ P,-) 2g T fB7 f ■ (30) 



For the calculation with a 1 resonance B\ as intermediate state in Fig.|3|b, following ref.pT 
we need the coupling between the p, the B and this resonance. This coupling has the 
following form 

Cbb iP = g s m Bl BB>?p, + [B{d»B v x )d vPll - (d v B)(d^) Ptl ) , (31) 
which leads to the following matrix element for the transition B° — > B^ p + \ 

< B^p + \B° >= g.m Bl e* Bl e; - -^-[e* Bl ■ (p B + p p ))(e* p ■ p Bl ). (32) 

TTl Bl 

With these definitions for the two couplings proportional to g s (s-wave transition) and g<i 
(d-wave transition) the weak transition matrix element becomes 



< n-p + \H eff \B° >= ^V ub V: d a 2 U(e* ■ Pw - 



m B — m p 
9sf Bl o 9df Bl 



m 2 Bl 



(33) 



In analogy to ( p9[) we relate the strong couplings in (|33|) to the transition from factors ^i(O) 
and ^2(0) appearing in the current matrix elements of the axial- vector current as defined in 

9sf Bl = (m B + m p )Af^ p (0), 

9df Bl = m ^ A B ^(0). (34) 
[m B + nip) 

Then (|33| ) can be written in terms of A\ and A2 or Aq: 

< n-p + \n eff \B° > = -Lv uh V^a 2 U{e* ■ p n -) \{m B + m p )A^ p (0) - (m B - m p )A^ p (0) 
= |Wa2/,(e*'A-)2mX" P (0). (35) 
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In (p3|) the pole at q 2 = m\ = m 2 Bi has canceled, as to be expected. The final result (j35l) 
differs from (|28|) only due to the cancelation of the pole, so that Aq^ p is to be evaluated at 
q 2 = instead of q 2 = m 2 in (|28|) . It is clear that this change makes only a small difference 
since A^ p (ml) ~ <^ p (0). 

The relation between the two different calculations of the pole diagram in Fig.|3|b is 
quite analogous to the results obtained for B° — > 7r + 7r~ in the previous section. We obtain 
analogous results for the decays B° — > D* + D~, D* + D~ , D* + n~ and D* + K~. Of course, for 
some of these channels the results ( p8|) and (^) differ by a larger amount, since in (^) the 
form factor Aq^ d must be calculated for a larger g 2 instead of q 2 = in (|35|). 

Similarly to the decay channel B° — > D + tt~ we have also contributions of the annihilation 
diagram for the channels B° — > D + p~ and B° —>■ D* + tt~ . If we evaluate them with the usual 
BSW approach or in the pole model with intermediate 0~ poles we obtain a contribution 
proportional to AQ^ p (m 2 B ) for B° — > D + p~ and Aq* ^ (m 2 B ) for the channel B° — > D* + 7r _ . 
These form factors, which must be evaluated at q 2 = m 2 B , are suppressed. Therefore we 
shall neglect the annihilation contributions as we have done also for the B° — > D + n~ decay. 



Only when we follow ref.[|TTJ, where the pole term is canceled by the same factor in the 
numerator, we arrive at a somewhat larger result which can be expressed again by ^^"^(O) 
and ^^'"^(O), respectively. For example, for the channel B° — > D + p~ the final result is: 

< D + p-\H eff \B° >= -^V cb V: dai f B (e* ■ p D+ )2mX" + "(0). (36) 

To derive this result we used the strong coupling Lagrangian for the coupling of the 1 + 
resonance D® to D + p~ and then introduced A^ p and A®^ p with (|34|) . A similar formula 
as (|36|) can be derived for the annihilation term of the decay B° — > D* + n~ . 

In the same way we calculate the class I decays of the charged B~ . These are the channels 
B~ — > D°D*~ and D°D*~ , which have amplitudes proportional to Ff^ D , and the channels 
B~ — > D*°D~ and D*°D~ with amplitudes proportional to Aq^ d * . 

The class II decays with amplitudes proportional to ci2 are calculated in the same way. 
The amplitudes are either proportional to Fi or to A Q . The decays B° — > 7r°J/^, K°J/ip 
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Table 8: B — >• PV class I decay amplitudes(10 7 ) and branching ratios in two versions of pole 
models and experimental measured branching 







Pole 


model 






BSW 


model 






Experimental 


Channel 


Amplitude 


BR 




Amplitude 


BR 




BR 


B° - 


-> TT + p 


0.033 a 2 


2.84 x 


10' 


-o 


0.033 a 2 


2.84 x 


io- 


c 

-o 


< 8.8 x 10~ 5 


B° - 


-> D+p- 


0.73a 2 -0.44ai 


6.60 x 


10' 


q 
-o 


0.73 a 2 


9.21 x 


io- 


q 

-o 


(7.8 ± 1.4) 10" 3 


B°- 


-> D+K*~ 


0.19 a 2 


4.62 x 


io- 


A 

-4 


0.19 a 2 


4.62 x 


io- 


A 
-4 




B° - 


-> D+D*- 


0.53 a 2 


3.72 x 


10" 


A 

-4 


0.53 a 2 


3.72 x 


io- 


A 

-4 


< 1.8 x 10~ 3 


B° - 




2.88 a 2 


8.90 x 


10" 


-3 


2.88 a 2 


8.90 x 


io- 


-3 


< 2.0 x 10~ 2 


B° - 


-> 7r _ /9 + 


0.016 a 2 


7.01 x 


10" 


-D 


0.016 a 2 


7.02 x 


10" 


-o 


< 8.8 x 10" 5 


B°- 




1.02a 2 -1.15ai 


1.37 x 


10" 


q 
-o 


1.02 a 2 


2.71 x 


io- 


q 

-o 


( 2.81 ± 0.24)10~ 3 


B°- 




0.28 a 2 


1.98 x 


10" 


A 
-4 


0.28 a 2 


2.01 x 


io- 


A 
-4 




B°- 


-> D~D* + 


0.43 a 2 


2.45 x 


10" 


-4 


0.48 a 2 


2.95 x 


io- 


-4 


< 1.8 x 10~ 3 


B°- 


-> DjD* + 


2.20 a 2 


5.75 x 


10" 


-3 


2.44 a 2 


7.09 x 


io- 


-3 


( 1.2±0.6)10 -2 


B~ - 


-> D°D*- 


0.53 a 2 


3.88 x 


10" 


-4 


0.53 a 2 


3.88 x 


io- 


-4 




B~ - 


-> D°D*~ 


2.88 a 2 


9.29 x 


10" 


-3 


2.88 a 2 


9.29 x 


io- 


-3 


(1.2 ± 1.0)10 -2 


B~ - 


-> D*°D- 


0.43 a 2 


2.56 x 


10" 


-4 


0.475 a 2 


3.07 x 


io- 


-4 




B~ - 


-> D*°D- 


2.20 a 2 


6.01 x 


io- 


-3 


2.44 a 2 


7.37 x 


io- 


-3 


(10 ± 7)10~ 3 



and B — > 7r J/ip, K J ftp belong to the first group. The second group consists of the 

decays B° — > p rj c , K*°rj c and B~ — > p~r/ c , K*~r] c . Other class II decays are the channels 
^ D ^ D o p o and go ^ p o n o_ 

The class III decays with amplitudes coming from both operators 0\ and 2 are the fol- 
lowing: B~ — > D°p~, D°K*~ and 5~ — > D*°tt~ , D*°K~ originating from 6 — > c transitions. 
The channels of class III with 6 — > it transitions are I? - — > p°7r~ and p~7r°. 

We have calculated all these decays again in the two versions, (i) pole model without 
canceling of the pole of the intermediate state, which gives the same results as the usual 



BSW calculation and (ii) pole-model in the form as introduced in ref. ||ll|| . This model has 
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a larger annihilation contribution and in the direct contribution, where the 0~ intermediate 
state is replaced by the 1 + state, the form factor is evaluated at q 2 = instead of q 2 = m 2 
or q 2 = m 2 D as it appears in version (i). 

The results for the class I decays are collected in Table 8 and for the class II and class 
III decays in Table 9. In these two tables, the amplitudes are defined without the factor 
of (e ■ p n ) in (|30"D. In Table 8, four of the decay channels can be compared to experimental 
data. In all four channels the agreement is rather good. Unfortunately the experimental 
errors for the £>(*)£>« final states are rather large. For B° — > n D* + the error is smaller 
and only the prediction of the BSW pole model agrees with the experimental number. The 
prediction of the pole model with canceling residues is a factor of two smaller than the 
experimental branching ratio. This is due to the additional annihilation contribution in this 
model which is too large when compared with the experimental branching ratio. Since in 
this model the off-shell extrapolation of the form factor to the q 2 = m 2 B is neglected the 
discrepancy with the data is not unexpected. Actually the disagreement of this pole model 
in the case of B° — > n~D* + is a strong argument against this version of pole model. Of 
course, this depends on the values assumed for ai and i.e. our assumption N = 2 also 
in this pole model. If a\ would be decreased as it happens for example, with N = 3, i.e. in 
pure factorization with colour-octet terms omitted, the annihilation term can be suppressed. 
However, as we shall see later, this leads to strong disagreement in other channels, where 
the two pole model versions have identical results. From this comparison it is clear that the 
extra form factor in the annihilation diagram must be taken into account also in this pole 
model. This makes these contributions negligible also. This decay gives us a good value for 
Aq~* d *(0). We notice also that the prediction for B° — > D~D* + agrees in the BSW-pole- 
model slightly better with the experimental number than in the other pole model. In this 
case the different results in the two models come only from the usual W emission graph. 

As already mentioned above, the two pole models give identical results for the decay 
channels B° — > D + K*~ , D + D*~ , D + D*~ and 7i + p~ . Differences occur for the decays B° — > 
D-D* + and D;D* + . 
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Table 9: B — > PV class II and III decay amplitudes(10 7 ) and branching ratios in two 
versions of pole models and experimental measured branching ratios or upper limits (90% 



c.l.) m, m, n 







Pole model 




BSW model 




Experimental 


Channel 


Amplitude 


BR 




Amplitude 


BR 




BR 


B°- 


-> vr° J/ip 


0.58 oi 


1.09 x 10" 


-5 


0.58 a t 


1.09 x 10" 


-5 


< 6.9 x 10~ 3 


B°- 


-> if j/v> 


4.07 ai 


4.89 x 10" 


-4 


4.07 ax 


4.89 x 10" 


-4 


(8.5t{J)10- 4 


B- - 


-> ir~J/ip 


0.83 ai 


2.25 x 10" 


-5 


0.83 ax 


2.25 x 10" 


-5 


< 4.4 x 10" 5 


B- - 


-> if- J/V> 


4.07 ai 


5.09 x 10" 


-4 


4.07 ax 


5.09 x 10" 


-4 


(10.2 ± 1.1)10" 4 


B° - 


-► 


0.08 ai 


3.34 x 10" 


-6 


0.12 ax 


7.20 x 10" 


-6 


_ 


B°- 


- if *V 


0.69 oi 


1.58 x 10" 


-4 


0.99 ax 


3.30 x 10" 


-4 


— 


Br - 




0.12 oi 


6.94 x 10" 


-6 


0.17 ax 


1.49 x 10" 


-5 




B- - 


if *-r/ c 


0.69 ai 


1.64 x 10" 


-4 


0.99 ax 


3.42 x 10" 


-4 




B°- 




0.024 ai 


1.03 x 10" 


-6 


0.024 ax 


1.03 x 10" 


-6 


< 2.4 x 10" 5 


B°- 


-» vr°i)* 


0.84 ai 


1.20 x 10" 


-4 


0.84 ax 


1.20 x 10" 


-4 


< 9.7 x 10~ 4 


B°- 




0.23 oi 


6.22 x 10" 


-5 


0.27 ax 


8.01 x 10" 


-5 


< 5.5 x 10~ 4 


B~ - 




0.02o 2 +0.01ai 


1.87 x 10" 


-5 


0.01a 2 +0.02ai 


1.87 x 10" 


-5 


< 7.7 x 10" 5 


B~ - 




0.01a 2 +0.02ai 


8.39 x 10" 


-6 


0.01a 2 +0.02ai 


8.39 x 10" 


-6 


< 4.3 x 10~ 5 


B~ - 




0.73a 2 +0.33ai 


1.19 x 10" 


-2 


0.73a 2 +0.38ai 


1.19 x 10" 


-2 


(1.34±0.18)10~ 2 


B- - 




1.02a 2 +1.19ai 


4.74 x 10" 


-3 


1.02a 2 +1.19ai 


4.74 x 10" 


-3 


(4.34±0.51)10~ 3 


B- - 


-> if-£>*° 


0.28a 2 +0.32ai 


3.42 x 10" 


-4 


0.28a 2 +0.32ai 


3.45 x 10" 


-4 




B- - 


-> if*-D° 


0.19a 2 +0.10ai 


6.15 x 10" 


-4 


0.19a 2 +0.11ai 


6.35 x 10" 


-4 
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In Table 9, our results for the branching ratios of the decays B° — > K°J/ip and B — > 
K~J/ip are somewhat smaller than the measured branching ratios. The factor r needed to 
achieve agreement is r = 1.7 ± 0.3 for B° — > K°J/%p and r = 2.0 ± 0.2 for B~ — > K~J/ifj, 
respectively. Since the Wilson coefficient a\ is already constrained by the class III decays 
the disagreement means either a (30 - 40)% larger value for Ff^ K (0) or a breakdown of 
the factorization assumption ||17|| . It is clear that other values for a±, which for example for 
iV = 3 in the pure factorization approach with colour octet terms neglected are obtained, 
would give very bad results for both pole models. For some of the other decays of class II 
upper limits for the branching ratios exist which are all obeyed by the theoretical results. 

The prominent class III decays are B~ — > p~D° and B~ — > tt~D*°. The predictions agree 
quite well with the experimental branching ratios. It is clear that this agreement, given the 
matrix elements of BSW, is only possible, if a\ is positive, as was already found in ref. |L6|. 



This statement is again independent of the form of pole model applied in the calculation. 
Concerning matrix elements these two decays depend on Ff^ D , Ff^ n , Aq^ p and Aq^ d * . 

Next we consider the B — > VV decays. For these decays the two pole-model versions 
give identical results which are the same as in the quark diagram approach. So, we need not 
write down all the formulas. We do this only for one channel B° — ^ p + p~ to fix our notation. 
The axial part of the current matrix element was already written in (|26|). In addition we 
need also the matrix element of the vector current which we write as usual 

< p + \u lfl b\B° >= ^ U *PBP" P { ^fl p y (37) 

With this, (|9|) and (|26"D we obtain for the weak transition matrix element 

< p + p-\H eff \B > = ^V ub V: d (-if p m p ) {e*(p + )e*(p-)(m B + m p )A^ p (m 2 p ) 

~ (e*(p + )- PB )(e*(p-)- PB ) 2A ^ P{ml) 



(fflg + m 



pi 



2V B ~* p (m 2 ) ) 

- *w^(p-K*(P + K^ +7 — -Mf[. (38) 

p (m B + m p ) J 

So, the transition matrix element is a linear combination of the form factor Ai, A-i and V 
which determine the transition in the s, d and p waves respectively [[|7] . 
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Table 10: B — ► VV decay amplitudes(10 8 ) and branching ratios in pole model and experi- 
mental measured branching ratios or upper limits (90% C.L.) fl6 , |2l 









Amplitudes 




model 




Experimental 


Channel 


v(ht 8 ) 


Al(10~ 8 ) 


A2(10~ 8 ) 


BR 




BR 


B°- 


-* P + P 


1.94 a 2 


0.84 a 2 


1.67 a 2 


2.40 x 10" 


-5 


< 2.2 x 10~ 3 


B°- 


-» D*+p- 


54.1 a 2 


25.0 a 2 


52.6 a 2 


8.70 x 10" 


-3 


(7.3 ± 1.5)10- 3 


B°- 


-» D*+K*~ 


14.3 a 2 


6.59 a 2 


13.9 a 2 


4.68 x 10- 


-4 


— 


B°- 


-> D* + D*- 


39.8 a 2 


18.1 a 2 


38.2 a 2 


9.28 x 10" 


-4 


(5.3^|)10- 4 


B°- 


-> D*+D*- 


215 a 2 


97.8 a 2 


206 a 2 


2.46 x 10" 


-2 


(1.9±1.2)10 -2 


B~ - 


-> D*°D*- 


39.8 a 2 


18.1 a 2 


38.2 a 2 


8.87 x 10- 


-4 




B- - 


-> D*°D*- 


215 a 2 


97.8 a 2 


206 a 2 


2.56 x 10- 


-2 


(2.3 ± 1.4)10~ 2 


B°- 


-> p° J/V> 


34.8 oi 


14.9 ai 


29.7 a x 


5.16 x 10- 


-5 




B°- 


-» K*°J/ip 


245 ai 


101 oi 


204 ai 


2.03 x 10- 


-3 


(1.32 ±0.24)10~ 3 


B- - 


^p-J/i; 


49.3 oi 


21.0 ai 


42.0 ai 


1.07 x 10- 


-4 




Br - 


-> K*~J/^j 


245 ai 


101 oi 


204 oi 


2.11 x 10- 


-3 


(1.41 ±0.33)10~ 3 


B°- 


^p°p° 


0.97 ai 


0.42 oi 


0.84 oi 


7.74 x 10- 


-7 


< 2.8 x 10~ 4 


B°- 


-> p°D* 


50.3 oi 


21.6 ai 


43.1 oi 


2.40 x 10" 


-4 


< 1.17 x 10~ 3 


B~ - 




1.37 (a 2 + ax) 


0.59 (02 + oi) 


1.18 (a 2 +ai) 


1.95 x 10- 


-5 


< 1.0 x 10~ 3 


B~ - 




54.1a 2 +71.1ai 


24.9a 2 +30.5ai 


52.6o 2 +61.0ai 


1.37 x 10- 


-2 


(1.55 ±0.31)10~ 2 


B~ - 


-> K*-D*° 


14.3a 2 +18.4ai 


6.59a 2 +7.95ai 


13.9a2+16.1ai 


7.50 x 10- 


-4 
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We have calculated the branching ratios for various decays using the form factor values 
given in Table 3 together with the decay constants in Table 2 and the pole masses from 
Table 1. The results are given in Table 10 and are compared to the experimental branching 
ratios (last column in Table 10). The decay channels in Table 10 are ordered according to 
class I, II and III as in the previous tables. The calculated branching ratios agree well with the 
experimental data, in particular for the decays B° — > D* + p~ , D* + D*~ and B~ — > D*°D*~ 
which are class I decays and B~ — > p~D*° which is a class III decay. The class II decays 
B° — > K*°J/tp and B~ — > K*~J/tp come out somewhat larger than the measured branching 
ratios. This could be changed by adjusting the B — > K* form factor by a small amount or 



may be a sign of the breakdown of the factorization assumption for class II decays (JT7 . 

In Table 10 we have collected also the amplitudes proportional to the form factors V B ^ P , 
A B ~^ P and A B ^ f '. The numbers in the first, second and third column of Table 10 are defined in 
such a way that the polarization dependent factors [— e^a^* (p~)e u * (p + )p B Pp+ / (rriB +m P ) 2 ] 
in the term proportional to V in fl3"8|), [e*(p + )e*(p~)] in the term proportional to Ai in ( |3"8"D 
and [(e*(p + ) • pb)(c*(p~) ■ Pb)/{itib + rn P ) 2 ] in the term proportional to A 2 in (37) are not 
included. Except for these factors the amplitudes are of the same order of magnitude as to 
be expected from the values of the form factors in Table 3. The results in Table 8, 9 and 10 
depend on the parameters A B ^ P , A B ~^ P ', V B ^ P and Af^ D * , A B ^ D * , v B ^ D * , respectively. 
The BSW values can be compared to results from other sources, i.e. either other quark 
models, QCD sum rules or lattice gauge theory. Values for A B ^ P , A B ~^ P and V B ^ P , reported 
in the literature, are collected in Table 11. Except for a few cases, the form factors are similar 
to the BSW values. From lattice gauge theory calculation the following values for Af~ >p (0) 
are reported: 0.22 ±0.05 |3§, 0.24 ±0.12 |3§ and 0.27^ 0, for a recent review see |!TJ. 



A similar collection for Af^ D \ A B ^ D \ V B ~* D * is found in Table 12. A B ^ P and A B ~* D * 
can be calculated with the help of ( p7|) from Ai and A 2 , respectively, given in Table 11 and 
Table 12. 

It is clear that the good agreement of the calculated partial decay rates with the exper- 
imental data in the case of the B — > PV and B — > VV decays serves only as a consistency 
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check of the form factors A\, A2 and V for B — > p and B — > D*. In order to obtain 
complete information on these form factors many more measurements are needed than just 
the branching ratios. Additional information can come from decay angular distributions of 
subsequent decays of the vector mesons. This has been studied in previous work W7\, and 



compared to experimental data in [16 



Table 11: Form factors Af^ p , A B ^ f ', V B ^ P at q 2 = in different models. 



Reference method 



m 



[0 
n 




c 
c 
c 
c 
a 
a 
a 
b 



Af^ p (0) 



0.24 ±0.04 

0.5 ±0.1 
0.35 ±0.16 
0.27 ±0.05 

0.283 

0.27 

0.26 

0.21 



A B - p 



(0) 



0.4 ±0.2 
0.42 ±0.12 
0.28 ±0.05 

0.283 

0.30 

0.24 

0.20 



V B ^P{0) 



0.28 ±0.08 

0.6 ±0.2 
0.47 ±0.14 
0.35 ±0.07 

0.329 

0.29 

0.35 

1.04 



We are now in the position to calculate other strong coupling constants of the B mesons, 
which enter into the pole-model formulation of non-leptonic weak decays, on the basis of the 
information following from the form factors A\, A2 and V for B — > p, B — > K* and B — > D*. 
As a reference we employ the BSW results for these form factors as given in Table 3. It 
is clear that the calculation can be repeated for any other choice as collected in Table 11 
and 12. Let us start with gsBp- This follows from (|29|) with the result gsBp — 1.54. In the 
same way we obtain gBB a K* = 1-81 and (}BB a D* = 3.48. We observe a large splitting of the 
coupling between two pseudoscalar B mesons (BB, BB S , BB C ) and the vector mesons p, K* 



and D*, respectively. This is mainly due to the mass factor of the vector meson in (|29). 
The other coupling constants which follow from Ai(0) and ^(0) with the relations (34) 
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Table 12: Form factors A B ^ D * , A B ^ D * , V B ^ D * at q 2 = in different models. 



Reference method 


Af^ D *(0) 


A^*(0) 


V B ^ D *(0) 


[ 

1 

I 


31 

i 

38 
39 


c 
a 
a 

f 


0.46 ±0.02 
0.651 
0.62 
0.622 


0.53 ±0.09 
0.686 
0.61 
0.778 


0.58 ±0.03 
0.705 
0.63 
0.747 



are the coupling constants of axial vector mesons Bi, B s i and B c \ with Bp, BK* and BD*, 
respectively. From V(0) we obtain the coupling constants of the vector mesons B*, B* and 
B* with Bp, BK* and BD*. The formula for this coupling in the case of Bp in analogous 
to the second formula in ( |34"| ) 

g = f B * V B ^(0). (39) 

The result of the coupling of the axial vector mesons is: 

g s (BB lP ) = 10.7, g d (BB lP )=8A, 
g s (BB sl K*) = 11.2, g d (BB sl K*) = 10.3, 
g s (BB cl D*) = 14.8, g d {BB cl D*) = 13.3. 

The corresponding vector meson coupling constants are: 

g(BB*p) = 9.64, g(BB*K*) = 9.67, g(BB*D*) = 12.1. 

We observe that all coupling constants have similar values, of the order of 10. The splitting of 
the couplings to Bp, BK* and BD* is somewhat larger than obtained for the corresponding 
VPP couplings to Bit, BK and BD in section 2. The coupling constant g(BB*p) can be 
compared with a QCD sum rule computation of this coupling |54 |. Using the values for f B 



and fs* as in our Table 2 these authors obtain g(BB*p) ~ 12 in reasonable agreement with 
our values above. 
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4 Summary and Conclusions 

We have calculated the branching ratios for the decays of neutral and charged B mesons 
into PP, PV and VV mesons using two versions of pole-dominance model in addition with 
a factorization assumption. The first pole model gives the same result as the extensively 
studied factorization model in terms of current matrix elements, if these current matrix 



elements are approximated by single poles. The second pole model due to Bedaque et al. |TT 
employs a different off-mass- shell extrapolation of the residues of the single-pole dominance 
approximation which leads, besides other differences, to an enhancement of the annihilation 
diagram contributions. From the derivation we see that the two versions of pole models 



are related. The pole model of ref.[|ll] involves a different way of extrapolation away from 
the pole position which has a particular large effect for some decays where an annihilation 
contribution is present. In the particular decay B° — > tt~D* + the model of ref. disagrees 



with the experimental branching ratios due to a large annihilation contribution. In the cases 
where annihilation diagrams do not contribute the final results of the two models do not 
differ very much. For the cases B — > VV and most of the B — > PV decays considered in 
this work they even give identical results. 

The coupling constants in the pole model are related to the usual form factors of the 
current matrix elements at q 2 = 0. For an overview we have compared the form factors Ff^ 71 
and Fi^ D and the form factors for the transitions B — > p and B — > D* calculated with 
different methods and found reasonable agreement. Using information for the pseudoscalar 
and vector meson decay constants we deduced the various strong couplings of the vector 
mesons B*, B* and B* with the Bn, BK and BD system, respectively. We found rather 
similar values for these three couplings showing that the splitting of the form factors Ff^ 71 , 
Ff~* K and Ff^ D is too a large extent related to the splitting of fs* and /#». The 
same pattern evolves for the couplings of the axial- vector mesons B\ and the vector mesons 
B* to the Bp, BK* and BD* system. It would be interesting to know whether the relations 
for the strong coupling constants found in this work could be obtained from QCD sum rule 
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or lattice calculations. 
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